Energy levels in polarization superlattices: a comparison of continuum strain models 
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A theoretical model for the energy levels in polarization superlattices is presented. The model 
includes the eflFect of strain on the local polarization-induced electric fields and the subsequent effect 
on the energy levels. Two continuum strain models are contrasted. One is the standard strain 
model derived from Hooke's law that is typically used to calculate energy levels in polarization 
superlattices and quantum wells. The other is a fully-coupled strain model derived from the ther- 
modynamic equation of state for piezoelectric materials. The latter is more complete and applicable 
to strongly piezoelectric materials where corrections to the standard model are significant. The 
underlying theory has been applied to AlGaN/GaN superlattices and quantum wells. It is found 
that the fully-coupled strain model yields very different electric fields from the standard model. The 
calculated intersubband transition energies are shifted by approximately 5-19 meV, depending 
on the structure. Thus from a device standpoint, the effect of applying the fully-coupled model 
produces a very measurable shift in the peak wavelength. This result has implications for the design 
of AlGaN/GaN optical switches. 

PACS numbers: 73.21.Ac, 71.20.Nr, 85.30.Tv, 73.20.At, 73.61.Ey, 77.65.Ly 
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I. INTRODUCTION 

Intersubband optical transitions (ISBT) in Al- 
GaN/GaN superlattices (SLs) and multiple quantum 
wells (MQWs) are being exploited for use in near- and 
mid-infrared lasers and ultra-fast all-optical switches in 
the 1.5-3 /zm wavelength range. liS*2*^i^*S4i A key de- 
sign issue related to ISBT-based device concepts is the 
calculation of the electron energy levels of these struc- 
tures so that the peak wavelengths can be estimated be- 
fore growth and fabrication. In the present work, we 
investigate the role of strain and polarization on the sub- 
band structure of SLs in the wurtzite crystal structure. 
Through a theoretical examination of fully-coupled and 
semi-coupled electromechanical treatments we show the 
importance of using a fully-coupled model for predicting 
the energy levels of SLs in strongly piezoelectric material 
systems. The model is then used to predict the energy 
levels of ten actual structures in order to compare our 
calculated ISBTs to previously measured spectra. 

Unlike zincblende semiconductor SLs, a number of is- 
sues arise in wurtzite SLs that complicate the task of cal- 
culating the energy levels. AlGaN in the wurtzite phase 
has a large spontaneous polarization moment along the 
[0001] axis. In addition, SLs grown on a SiC or sapphire 
substrate are pseudomorphic and the large in-plane bi- 
axial strains induce a piezoelectric polarization moment 
oriented along the c-axis with the direction depending on 
whether the strain is tensile (SiC) or compressive (AI2O3). 
The discontinuity of the polarization moments effectively 



represents fixed sheet charges at the interfaces of the SL. 
In general, each AlGaN on GaN interface in the direction 
[0001] will have a positive space charge and each GaN 
on AlGaN interface a negative space charge. Thus unlike 
zincblende SLs which are flat band unless doped, a calcu- 
lation of the electronic eigenvalues using the Schrodinger 
equation must be preceded by a calculation of the elec- 
trostatic potential, representing the Hartree term in the 
Schrodinger equation, using the Poisson equation. 

There is a further complication that has been pre- 
viously ignored in calculations of the energies in Al- 
GaN/GaN SLs. This involves the incorporation of strain 
into the electric field and eigenvalue calculations. To 
date, the strain model for AlGaN/GaN SLs has been bor- 
rowed from the zincblende realm;^ without additional con- 
sideration given to its validity for strongly piezoelectric 
materials. Although piezoelectric, zincblende materials 
have comparatively small piezoelectric tensor elements 
so that the thermodynamic equation of state is reduced 
to the standard Hooke's law with little or no error. From 
this relation, the strain tensor for zincblende SLs can be 
readily worked out with good accuracy. On the other 
hand, wurtzite materials have large piezoelectric coeffi- 
cients indicating strong coupling between the strain and 
electric fields. In the present case of group Ill-nitride ma- 
terials, we will show that treating the mechanical strain 
as separate from the electronic properties is no longer a 
sound methodology. The result is that a linear stress- 
strain model (Hooke's law) is no longer valid, and the 
fully-coupled thermodynamic equation of state must be 
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invoked to obtain the strain and electric fields simulta- 
neously. 

The coupling described in the present work is some- 
what analogous to the electromechanical coupling in sur- 
face acoustic wave (SAW) devices using AIN and GaN 
thin films pSiifiiiiiiSiii The strength of the interaction be- 
tween the electronic and mechanical properties in SAW 
devices is determined by the electromechanical coupling 
coefficient^ a quantity that measures the interaction be- 
tween the acoustic and electromagnetic waves in piezo- 
electric materials. In contrast, the electromechanical 
coupling described herein deals with the interaction be- 
tween the static electric and strain fields. Although the 
mathematical treatments of the two cases are very dif- 
ferent, both types of couplings originate from the same 
thermodynamic equation of state. In other work, the 
fully-coupled theory has predicted deviations in the static 
strain fields present in AlGaN/GaN heterostructure field- 
effect transistors"'^^ and the idealized case of free-standing 
(as opposed to substrate-conforming) superlatticesfli 

In previous modeling of piezoelectric SLs, the me- 
chanical and electronic properties are treated separately 
and sequentially: (i) first the in-plane strain is calcu- 
lated from the pseudomorphic boundary condition, (ii) 
Hooke's law is then invoked to obtain the longitudinal 
strain, and (Hi) the calculated strain tensor is subse- 
quently used as an input to the Poisson and Schrodinger 
equations. The strain is never recalculated to reflect the 
presence of static electric fields in the constituent lay- 
ers of the SL. In this paper, we compare the standard 
approach with a more rigorous continuum elastic theory 
applicable to piezoelectric materials. Using the proposed 
formalism, we apply the fully-coupled equation of state 
for piezoelectric materials to obtain simultaneously the 
strain and electronic properties of AlGaN/GaN SLs. 

First we treat the case of undoped SLs and show that 
closed-form analytical expressions can be obtained for 
both the strain and electric fields, following which the 
eigenstates can be calculated using the Schrodinger equa- 
tion. It will be shown, using specific examples, that 
the calculated strain and electric fields differ substan- 
tially from those obtained using the standard (uncou- 
pled) strain model. Depending on the Al fraction and 
the geometry of the SL, the longitudinal strain calcu- 
lated from the standard model may be in error by as 
much as 40% relative to the fully-coupled model. Fur- 
ther, it will be shown that the calculated ISBT energy 
may differ from that of the standard model by as much 
as 16meV, depending on the SL geometry. 

Second, we treat the more useful case from a device 
standpoint of SLs Si-doped in the well. (The doping 
can be tailored to populate the lowest conduction sub- 
band to facilitate optical transitions.) In this case, it 
is not possible to obtain closed-form analytical expres- 
sions for the strain and electric fields. Instead, we use a 
Schrodinger-Poisson solver in conjunction with the fully- 
coupled equation of state. The peak wavelength is calcu- 
lated for a number of structures and the results compared 



with published experimental data. Once again it will be 
shown that the standard and fully-coupled models yield 
significant differences in the ISBT energy. 

This paper is organized as follows: In Sec.^the contin- 
uum strain model is described. In Sec. Ill Al the general 
equations for the fully-coupled strain model are obtained. 
In Sec. lII BI the strain tensor and electric field for a polar- 
ization SL are worked out. In Sec. lII Gl the calculation of 
the electron eigenstates is described. A fully-coupled nu- 
merical model is outlined in Sec. Ill DI Galculated results 
are presented in Sec. lIIII In Scc. lIII Al calculated results 
for the standard and fully-coupled cases are contrasted 
for a model undoped SL. In Sec. IIII BI both models are 
tested against published experimental data for a series of 
doped SLs. The results are summarized in Sec. II VI 

II. MODEL DESCRIPTION 

Ordinarily, calculating the strain or stress tensor for 
a generalized strain problem becomes a complicated nu- 
merical exercise involving minimizing the Helmholtz free 
energy within the problem domainii^ This approach suf- 
fices for most materials, but specifically not for strongly 
piezoelectric materials. The reason can be illustrated 
as follows: If we take a piezoelectric plate and apply 
an external stress to it, the plate will be geometrically 
deformed and, because of the piezoelectric effect, a po- 
larization moment will be induced, accompanied by an 
internal electric field. But in addition to the piezoelec- 
tric effect, there is also a converse piezoelectric effect. 
In our plate example, the induced electric field result- 
ing from the external stress will exert a counter force 
to resist deformation of the plate. In a self-consistent 
way, the crystal will reach its equilibrium state consonant 
with minimum stored energy. This effect is present in all 
non-centrosymmetric crystals, but is especially strong in 
certain hexagonal crystals. Consequently, the uncoupled 
strain model for zincblende SLs results in errors when ap- 
plied to wurtzite SLs if the converse piezoelectric effect 
is substantial, as in our case. 

The relevant energy functional for piezoelectric mate- 
rials is the electric enthalpy H given hyi^ 

H = U -EB, (1) 

where E and D are the electric field and electric displace- 
ment, respectively, and U is the total internal energy 
(strain -I- electrostatic) given by 

^ = ^Cijknijlki + ^EijE^Ej, (2) 

in which Cijki is the fourth-ranked clastic stiffness ten- 
sor, Eij is the tensor form of the electric permittivity, ^ij 
is the strain tensor, and the indices i, j, k, and I run 
over the Gartesian coordinates x, y, and z. Summation 
over repeated indices is implied throughout. Accompany- 
ing the energy functional is the constitutive relationship 
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for the electric displacement. For piezoelectric materials, 
this is given, with the spontaneous polarization included, 
by the expression 

A = e^jk'^jk + e^JEJ + i^', (3) 

in which e^fe is the piezoelectric coefficient tensor and 
is the spontaneous polarizationiSSi For wurtzite materials, 
only the z component of exists because of the sixfold 
rotational symmetry of the [0001] axis. The first term in 
Eq. (PI is the piezoelectric polarization. After substitu- 
tion into Eq. the final form of the electric enthalpy 
becomes 

^ = ^Cijknijlki - eijkEi'^jk - ^EijEiEj - EiP^. (4) 

A. Fully-coupled strain tensor for planar strain 

In principle, minimizing H within the problem domain 
gives the strain and electric fields for a generalized prob- 
lem. In practice, this often means having to set up com- 
plicated finite element calculations. Problems involving 
two- and three-dimensional geometric variations will be 
subjects of future numerical work and we instead focus 
here on the SL case where the issues of grid and mini- 
mization technique will not obscure the physics. The SL 
problem is a planar one-dimensional (1-D) strain problem 
with, at least nominally, no shear strains. Accordingly, 
we can begin from the linear piezoelectric equation of 
state 

= CijkHki — SkijEk, (5) 

obtained by differentiating Eq. I^J with respect to the 
strain tensor, where cr.y is the stress tensor. Expanding 
Eq. Q and using the Voigt notation^ for the third- and 
fourth-ranked tensors, the following stress-stain relation- 
ships are obtained, assuming the z axis to be the sixfold 
axis of rotation: 

c^xx ^IxxCii +7yyCi2 +72^^13 "Csi^^z, (6a) 



CTyy = ^xxCl2 + "fyyCll + JzzCl3 - e^lE^, (6b) 



In the absence of the electric field, these equations are 
recognized as the familiar tensor form of Hooke's law for 
hexagonal crystals. 

In conjunction with Eq. I^, we use the constitutive 
relations obtained by expanding Eq. (O: 

D, = P, + eE,, (7a) 



Dy^Py+eEy, (7b) 



Dz = Pz+eEz+P^, (7c) 

where the electric permittivity is taken to be isotropic, a 
reasonable approximation for AlGaN/GaN SLs, and the 
piezoelectric moments are given by 

Px = 2ei572;2, (8a) 



Py = 2ei57j/z, (8b) 

and 

Pz = esiijxx + jyy) + 6337^^. (8c) 

For simplicity, it is assumed that there are no shear 
strains, manifested by warping, within the structure. The 
boundary condition for a free surface, aiz — 0, can then 
be applied throughout the layers, instead of just at the 
surface. From Eq. l(Hc|) . this gives 

2Ci3 , 633 

Izz = T^lxx + -T^Ez, (9) 

L'33 L.33 

where "fyy = ^xx m the 1-D planar case and ^xx is as- 
sumed to be known from the pseudomorphic condition 
across the interfaces. There still remains the problem of 
finding the electric field which is the topic of the next 
section. 



{ixx + lyy)Cis. + lzzCs.3. - CzzEz, (6c) 



B. Poisson equation 



For an isolated piezoelectric plate under planar stress, 
the constitutive equations and the equations of state 
a^y — "ixy{Cii — C12), (6d) should be sufficient for obtaining the strain and elec- 

tric fields. For the SL, however, the continuity of the 
electric displacement must satisfied at the interface, and 
periodic boundary conditions must be imposed on the 
electrostatic potential (p, as well as the continuity of (p 
g^jjj across the interface. Additional complications will arise 

from doping, as this will give rise to space charges and 
ayz = 2jyzC4i — ei^Ey. (6f) free electrons. These requirements are all met by solving 



<Jxz = '^IxzCiA - ei5Ex, (6e) 
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the Poisson equation. From Gauss's law and Eqs. (f7c|l . 
(|5c)l . and lO, we obtain the 1-D Poisson equation 



d 



dz 



dz 



631 - 6337;— ) Ixx 
'-'33 



(10) 



where N'^ is the ionized donor concentration, n is the free 
electron concentration calculated from the Fermi energy 
and the wave functions, e is the electronic charge, and 



'=^33 
C33 



(11) 



It is evident from Eq. (|10|l that k serves as an effective 
electric permittivity in the fully-coupled case. Also, be- 
cause e§3/C33 > 0, the electromechanical coupHng results 
effectively in additional dielectric screening. 

Figure ^ shows the band edges for a period of the 
SL under consideration. In the following derivation, "a" 
refers to the AlGaN layer and "6" the GaN layer. Equa- 
tion (|10|l is solved to obtain subject to the continuity 
of (j) and the electric displacement Dz aX z ^ Wa- The 
latter is expressed by 



dz 



ps 



631 — 633 



^xx 



(12) 



It is assumed that there is no applied bias. Periodic 
boundary conditions then apply. This is accomplished 
by setting (j) = aX z — Q and z — {wa + wi,). Unless 
the free electron distribution can be realistically approx- 
imated by a ^ function, Eq. (|10|) has to be solved numer- 
ically in the most general case. To illustrate the concept 
of electromechanical coupling, we assume for the moment 
that the SL is nominally undoped and with no free elec- 
trons from traps or surface states. Later, we will present 
results for doped SLs using our fully-coupled numerical 
model. For a depleted SL, the general solution of Eq. 
((Tn|l is given by 

, , 2(e3iC33 - e33Ci3)7r.r. n^A 
= — z-\ — zH z + B, (13) 

K KG33 K 

where A and B are unknown constants. Thus there are 
four unknowns, two in each layer. All four constants are 
accounted for by the four boundary conditions discussed 
above. 

After obtaining the unknowns, the electric fields in the 
two layers are given by 



Wb 



(ps(b)_ps(a)) ^ 2wb{e-,Cf,-el,CI,)j-, 



WaKb + WbKa 



C^^{WaHb -I- WbHa) 



2«;b(e|3C?3-eb,C|'3)7b^ _ Wb 



C^-^iWaKb +WbKa) 



In the standard model, E'^ and i?^ are obtained by re- 
placing K by e in Eq. H14|l . using the appropriate sub- 
scripts for the two layers. It is seen, therefore, that the 



fully-coupled electric field is smaller than its standard 
counterpart. From Eq. lO, the longitudinal strain in the 
"a" layer is given by 



2C- 



13 



a 



IXX ' 



2^&^33 (^33^13 



-31^33)7^ 



CI3 [Wanb +WbKa) 



'^33^33i'^aKb + WbKa) 

and in the "6" layer by 



Cl^iWaKb +WbHa) 



-(15) 



I zz 



2C^ 

'-'33 



7*^ 

I XX 



2^0633(633(7- 



=31 



^33/ IXX 



C^g [WaKb + WbKa) 



2^«a6^3(e|3C: 



13 



liCf3)7 



w, 



^h^^Ssi'^o.'ib + WbKa) 



C^^iWaKb +WbKa) 



(16) 



We can compare these expressions for strain in the 
wurtzite system directly with the zincblende case where 
the spontaneous polarization terms vanish and the com- 
pliance tensor has fewer unique elements. There, a simi- 
lar (but somewhat less complicated) expression to those 
in Eqs. ifTBjl and lfTB|l is obtained for the longitudi- 
nal strain in a [lll]-oriented pseudomorphic layer. The 
zincblende [111] case was derived seperately by Bahder— 
using the method of Lagrange multipliers to minimize 
the free energy density, an alternate approach. In other 
work on lattice dynamics in undoped GaN/AlN SLs^S 
comparable electric field corrections to the strain along 
the growth direction are obtained with the main differ- 
ence being the use of the high-frequency dielectric per- 
mittivity e(cxD) as opposed to the present case of static 
screening, as in Eq. (fTT|l . 

The in-plane strains are calculated by assuming perfect 
in-plane atomic registry of the SL layers with the buffer 
layer. Applying this condition, the in-plane strains are 
given by 



Obfr 



and 



7"^ 

/ XX 



flbfr — a-b 



Gb 



(17) 



(18) 



where aa and ab are the relaxed c-plane lattice constants 
of the "a" and "&" layers, respectively, and Cbfr is the c- 
plane lattice constant of the buffer layer. The foregoing 
model also works for less than perfect registry: if the 
in-plane strains are known independently, they can still 
be substituted into the above equations to obtain the 
electric fields and longitudinal strains. As is well known, 
the standard model gives the longitudinal strains as 



a (std) 

/ zz 



2C: 



13 



C: 



Ix 



33 



and 



7 



b (std) 



^^^13 -,b 

^b '^a; 
'-'33 



(19) 



(20) 
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i.e. the first terms in Eqs. (|15|l and (|16|l . and, as a conse- 
quence, omits a great deal of physics under certain con- 
ditions. It will be seen shortly that the fully-coupled 
correction to the standard longitudinal strain is quite sig- 
nificant. 



C. Schrodinger equation 

Owing to the large band gaps of the consituents of 
the AlGaN/GaN SL, the electron eigenstates states can 
be described by a Hamiltonian in the T^^ basis without 
including any mixing from the Fg^ and Tiy hole states, 
incurring little error in the process. The resulting Hamil- 
tonian is the one-band Schrodinger equation 

2 dz \m* dz J 2m* 
{AEc -ec/j- e(/)^c)* + a^ilxx + lyy + Izz)"^ = -^(k)*, 

(21) 

where and ky are the electron wave vectors in the 
c-plane, ^ is the electron wave function, i?(k) is the to- 
tal electron energy, </> is the electrostatic potential dis- 
cussed in Sec. IIIBI and represents the Hartree part of 
the Coulomb interaction, (j)xc represents the exchange- 
correlation part of the Coulomb interaction, m* is the 
effective electron mass, Qc is the conduction band hydro- 
static deformation potential, AEc is the conduction band 
discontinuity before strain shown in Fig. ^ and ja has 
been defined previously. 

It should be noted that all of these quantities depend 
on z. For AE^ we assume that 60% of the band gap 
difference between the two materials appears in the con- 
duction band, with the caveat that the offset is not well 
known. One could legitimately use the conduction band 
offset as an adjustable parameter to try to fit published 
experimental data, but it has been kept fixed in the cal- 
culated results presented here. There is a net hydrostatic 
component of the strain obtained from the sum of the di- 
agonal elements of the strain tensor. This component will 
shift the band edge to higher or lower energy, depending 
on whether the hydrostatic component is compressive or 
tensile. 

If the SL is undopcd, the electric field is piecewise 
constant so that cj) — ~Fz in Eq. (|21|) in the respec- 
tive layers. Analytic solutions of the wave function then 
can be obtained using Airy functions. ■^^ A much more 
flexible approach, however, and the one adopted in the 
present work, is to use a discretized numerical technique, 
e.g. finite-differencing, that can also handle the more 
technologically interesting case of doped SLs. For SLs 
and MQWs, Bloch boundary conditions are enforced, i.e. 
vI'(O) = ^[wa + Wb) exp[— ifcz(wa + Wb)\^ whcrc kz is the 
crystal momentum corresponding to the periodicity of 
the layers along the growth axis. 



D. Fully-coupled numerical model 

For doped SLs in which the electrostatic potential is 
very non-linear, the Poisson and Schrodinger equations 
cannot be solved analytically in closed form. For this 
case, we use a fully-coupled numerical model. The cen- 
tral framework for this is a Schrodinger-Poisson solver. 
The electron states and the free electron distribution are 
calculated by solving Eq. H21|) on a finite-difference grid 
subject to the boundary conditions discussed above. If 
present, hole states are calculated using a 6x6 k-p Hamil- 
tonian. For the exchange-correlation potential, we use 
the parameterized expression of Hedin and Lundqvist,-4 
derived from density-functional theory within the local- 
density approximation. The charge-balance equation, 
which determines the position of the Fermi energy Ep 
in relation to the SL subbands, is solved by the Newton- 
Raphson method. Fcrmi-Dirac statistics are used for the 
probability of occupancy of the electron states. 

The model has been described in detail in Ref. 
and the band structure and strain parameters provided 
therein. Since then, a fully-coupled strain calculation 
has been added to the numerical model by solving the 
modified Poisson equation, i.e. using k instead of e as 
shown in Eq. and incorporating Eq. ^ into the 

self-consistent calculation. This means that the strain 
terms in Eq. (|21() are updated each time it is solved. In 
an uncoupled calculation, the strain terms would remain 
invariant throughout the self-consistent calculation. It 
has been shown'^^ that there is a bowing of the spon- 
taneous polarization as a function of x. This effect is 
included in the present model. 

III. RESULTS AND DISCUSSION 

First we show calculated longitudinal strains and elec- 
tric fields for a model undoped SL to illustrate the dif- 
ferences between the standard and fully-coupled strain 
models. We then show how these differences lead to dif- 
ferences in the calculated eigenstates of the SL. We then 
present calculated ISBT energies and peak wavelengths 
using both the standard and fully-coupled models for 
doped SLs and compare the results with published ex- 
perimental data. 

A. Undoped Superlattices 

We consider a model SL consisting of 20A A^Cai^^jN 
barriers and 60A GaN wells on a GaN buffer. Assum- 
ing the pseudomorphic condition to hold, the GaN layers 
will have no in-planc strain components, while the Al- 
GaN layers will have in-plane strains in accordance with 
Eq. lO. Using Eqs. lO, ijTsl), and (HBJ, the strain and 
electric fields are calculated for the fully-coupled model 
and compared with the standard results. Following es- 
tablished convention, a negative sign in the present cal- 
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culations indicates contraction and a positive sign exten- 
sion relative to the unstrained state. Figures|2Ja) and (b) 
show the longitudinal strain in the barrier and well layers, 
respectively, as a function of the barrier mole fraction for 
the fully-coupled and standard cases. Despite 7^2,=0 due 
to the lattice matching condition, a non-zero occurs 
due to electromechanical coupling, as predicted in Eq. 
© and again in Eq. and shown in Fig. ^h). This 

strain is a near linear function of the Al fraction and, in 
this example, is about —0.07% for x—1. 

The largest strains occur in the AlGaN layers due to 
the lattice mismatch and it is also here that a significant 
deviation between the standard and fully-coupled models 
is seen as shown in Figs.|2Ia) and (b). This deviation is 
shown in Figs.|3Ja) and (b). The error of the standard 
model relative to the fully-coupled model can be in excess 
of 35% AlN/GaN SLs, as seen in this example. It is 
even higher in structures with higher electric fields in the 
barrier. This occurs when Wb^Wa- For example, if we 
set Wa = W^ and u'6=60A, the error for a; = 1 is about 
45%. These deviations in the strain are quite significant 
and, as will be seen shortly, have an impact on the ISBT 
energies. 

Figure 0] shows the calculated electric fields in the Al- 
GaN and GaN layers for our model SL. From Eq. 
it is seen that the larger electric field occurs in the thin- 
ner layer. The electric field calculated from the fully- 
coupled model is smaller in magnitude than the standard 
electric field due to an effective screening caused by the 
electromechanical coupling. This screening increases at 
higher strains. The deviation between the standard and 
fully-coupled electric fields is about 7% for x = 1. 

Figure [S] shows the calculated ISBT energies between 
the first two electron subbands and the corresponding 
peak wavelengths for the model SL. The energies are cal- 
culated at kz{'Wa+Wb)—Tr, the location of the minimum 
energy separation between the first two subbands in the 
Brillouin zone. The present calculations show that there 
is little change in the energies between the zone center 
and zone boundary for a wide range of SLs. A num- 
ber of factors contribute to the relatively narrow mini- 
bandwidth. First, the band edge discontinuity AEc is 
quite large due to the large band gaps of the host mate- 
rials. Second, the effective electron mass is large, in this 
case, 0.2mo in GaN and 0.33mo in AIN. Third, the built- 
in electric field causes the electron wave function to be 
localized in the triangular notch close to the AlGaN/GaN 
interface (see Fig. All of these factors reduce the ex- 
ponential tail of the electron wave functions between ad- 
jacent wells, which, in turn, would appear as a dispersion 
in the mini Brillouin zone. For some SLs, however, par- 
ticularly those with thin wells, the wave functions will 
spread into the barrier layers, causing some dispersion in 
the Brillouin zone. 

The most significant feature of Fig.|5lis the discrepancy 
between the standard and fully-coupled models. For ex- 
ample, for X = 0.3, the fully-coupled transition energy is 
lower than the standard value by 3.7 meV. This differ- 



ence increases to 19.5 meV for x = 1. The latter result is 
especially significant, because high Al fractions are pre- 
ferred for optical switching technology due to the shorter 
peak wavelength. The difference in energies between the 
two models is large enough to be measurable by, for ex- 
ample, Fourier-transform infrared spectroscopy (FTIR). 
The wavelength for the standard model is shorter by 
about 4% relative to the coupled model for a; = 1. 

The reason for the red shift of the fully-coupled re- 
sults relative to the standard model can be understood 
by noting that the introduction of electromechanical cou- 
pling reduces the magnitude of the electric field (see Fig. 
0J| . The smaller electric field results in a conduction band 
profile closer to flat-band conditions (see Fig.QJ. In addi- 
tion, a lowering of the effective barrier height occurs from 
the reduction of 7^^ seen in Fig. [3 The more shallow tri- 
angular notch, together with a reduced barrier height, 
will cause the subbands to more closely spaced in energy. 
The calculated red shift of the ISBT is distinct from the 
Stark shift seen in interband transitions where the tran- 
sition energy shifts to higher energy as the electric field 
is reduced. 



B. Doped Superlattices 

For optical switching technology, it is necessary to n- 
dope the SL in order to populate the first electron sub- 
band to facilitate ISBTs. For such structures, we use the 
fully-coupled numerical model described previously. The 
model is tested against published optical data for various 
SL structures. TableQlshows the calculated ISBT energy 
between the first two subbands for ten SL samples taken 
from the literature. The standard and fully-coupled re- 
sults are contrasted. It is evident that two models give 
differing results. Also evident is the consistent red shift 
of the fully-coupled results compared to the standard re- 
sults for the reasons discussed in Sec. IIII Al The differ- 
ences depend on the layer thicknesses and doping of the 
samples, varying from 4.9 meV for sample A to 19.4 meV 
for sample D. These differences are significant enough to 
be measurable by standard techniques such as FTIR. 

Also shown in Tableware the experimentally-obtained 
peak wavelengths for the SLs. These are compared with 
the calculated wavelengths from the standard and fully- 
coupled models. Except for samples F and G, it is clear 
that the calculated wavelengths are in reasonably good 
agreement with the published data. The causes of the 
discrepancies for samples F and G are unclear at this 
point. It should be noted that we have not attempted to 
optimize the input parameters and have chosen instead 
to use a generic set of parameters'^ without fitting. The 
calculated results are very sensitive to all of the input 
parameters and also to the geometry and Al fraction. For 
instance, if the well thicknesses in samples F and G are 
increased by two monolayers and x reduced to 0.6, the 
wavelength can be fitted to within 5% using the fully- 
coupled model. More precise modeling of optical data 
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will be the subject of future work. For now, we simply 
wish to illustrate the importance of incorporating a fully- 
coupled strain model in the design of optical switches. 

The calculations in Table^were done for a temperature 
of 300K. At 77K, there is a blue shift of the transition 
energies due to the slight increase in Ai^c. The blue shift 
is largest for SLs with the thinnest wells wherein the 
subbands are pushed closer to band edge discontinuity 
and smallest for SLs with the thickest wells in which the 
first two subbands see less of the band edge discontinuity. 
For example, the shift is about 8.5 meV for sample A and 
about 0.15 meV for sample E. 

Figure shows the calculated conduction band edge 
and electron distribution for sample C in Table HI using 
the fully-coupled model. Also shown are the Fermi en- 
ergy and the first three electron subbands calculated at 
the Brillouin zone boundary. This profile was calculated 
at 300K. At 77K, there is no discernible change in the 
electron distribution function and the slope of the con- 
duction band edge. There are, however, shifts in the 
subbands of a few meV depending on the structure, as 
described earlier. As the calculation shows, the Fermi 
energy appears slightly above the first subband but well 
below the second subband, in spite of the high doping, 
ensuring that the first subband is populated by electrons 
and the second nearly empty in order to facilitate opti- 
cal absorption. This Fermi energy position is consistent 
with measured SL structures with transition energies cor- 
responding to Ei^E2 transitions. The calculated distri- 
bution and band edges, therefore, appear plausible. 

Figure [7| shows the electric field distribution for se- 
lected structures from Table using the fully-coupled 
model. The large electric fields in these structures are a 
consequence of the large polarization discontinuity across 
the interface. It is difficult to verify these fields directly. 
There is indirect evidence, however, that these fields are 
not unreasonable given the close fits of the ISBT wave- 
lengths with experimental data. Due to the heavy dop- 
ing, analytical expressions commonly used to estimate 
the electric fields would lead to errors, especially in the 



wells where the field is clearly non-linear. Even on the 
barrier side near the interface, there is an in increase 
in the magnitude of the field due to the penetration of 
the wave functions into the barrier. For such SLs, a nu- 
merical solution of the fully-coupled Poisson equation as 
describe here is essential. 



IV. SUMMARY AND CONCLUSIONS 

In summary, a fully-coupled model for the strain and 
the eigenstates of AlGaN/GaN polarization SLs has been 
presented. This model is compared with the standard 
strain model utilizing Hooke's law. Both the sponta- 
neous and piezoelectric polarizations are included, to- 
gether with free electrons and ionic space charges. It 
is seen that the strain and electronic properties of the 
material are linked through the fully-coupled thermody- 
namic equation of state for piezoelectric materials. Sepa- 
rating the mechanical and electronic aspects of the SL in 
any theoretical modeling of the properties of these struc- 
tures leads to errors in both the strain and the eigen- 
states of the system. For strongly coupled cases, such as 
AlGaN/GaN SLs, the corrections to the standard model 
can be significant. The ISBT energies calculated from 
the fully-coupled model show a measurable red shift com- 
pared to the corresponding energies calculated from the 
separable model. This result has consequences for the 
design of optical switches utilizing AlGaN/GaN SLs. 
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TABLE L Comparison of published intersubband data with 
our calculated results for Ali,Gai_2,N/GaN SLs. The calcu- 
lated energies in eV represent the separation between the first 
and second subbands at kz{wa+wt)—n. The peak wavelength 
A is in /xm. The experimental (exp) results are contrasted with 
calculated results from the standard (std) and fully-coupled 
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FIG. 1: Conduction and valence band edges of one period of 
a superlattice or MQW. Assuming that the substrate is to the 
right, the schematic depicts a cation-faced structure. Wa and 
Wb are the thicknesses of the AlGaN and GaN layers, respec- 
tively. The directions of the spontaneous and piezoelectric 
moments in the two layers are indicated, assuming the buffer 
is GaN. The dashed lines indicate the first three minibands 
and the dot-dashed line the Fermi energy. 
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FIG. 2: Calculated longitudinal strain in (a) the barrier layer 
and (b) the well layer for an undoped SL consisting of 20A 
Ala;Gai_iN and 60A GaN on a GaN buffer as a function of 
X. The fully coupled and standard results are shown. 
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FIG. 3: Difference between tlie fully coupled and uncoupled 
longitudinal strains in (a) the barrier layer and (b) the well 
layer for the SL of Fig.|21as a function of Al composition x in 
the barrier. 
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FIG. 4: Calculated electric field in (a) the barrier and (b) the 
well layers for the SL described in Fig. |5| as a function of Al 
composition x in the barrier. The fully coupled (solid lines) 
and standard (dashed lines) results are shown. 
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FIG. 5: Calculated ISBT energies (left y axis) and peak wave- 
lengths (right y axis) between the first and second subbands 
at kz{wa+Wb)=iT for the SL described in Fig.|21as a function 
of Al mole fraction x in the barrier. The fully coupled (solid 
line) and standard (dashed line) results are shown. 
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FIG. 6: Calculated conduction band edge (left y axis) and 
electron distribution (right y axis) in one SL period for sam- 
ple C in Table |I] using the fully-coupled numerical model. 
The first three electron subbands (dashed lines), calculated 
at kz{vLia+'Wb)=n , and the Fermi energy Ef (dot-dashed line) 
are shown. 




FIG. 7: Calculated electric field in one SL period for sam- 
ples A, C, and E in Table Q using the fully-coupled numerical 
model. The sign change in the electric field marks the position 
of the AlGaN/GaN interface. 



